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Abstract 

All eigenstates and eigenvalues are determined for the spin- 1 /2 XXZ chain 
H = 2JEi(SfS?+i + SiSf+i + ASfSf+i) for rin gs with up to N = 16 spins, 
for anisotropics A = 0, cos(0.3-7r), and 1. The dynamic spin pair correlations 
(Sj J +n (t)Sj l ) , (fi = x,z), the dynamic structure factors and the in- 

termediate structure factors I^(q,t) are calculated for arbitrary temperature 
T. It is found, that for all T, S z {q,uj) is mainly concentrated on the region 
M < £2(0), where £2(9) is the upper boundary of the two-spinon continuum, 
although excited states corresponding to a much broader frequency spectrum 
contribute. This is also true for the Haldane-Shastry model and the frustrated 
Heisenberg model. The intermediate structure factors I^(q,t) for A/0 show 
exponential decay for high T and large q. Within the accessible time range, 
the time-dependent spin correlation functions do not display the long-time 
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signatures of spin diffusion. 
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I. INTRODUCTION 



We discuss the one- dimensional S = 1/2 XX Z antiferromagnet, specified by the Hamil- 
tonian 

N 

Hxxz = 2Jj^ \SiSf +1 + SfS l - +1 + AS?S* +1 j , (1.1) 
i=i 

with coupling J > and anisotropy < A < 1. The model is defined on a ring of 
N lattice sites. The quantities studied in this paper are the dynamic spin pair correlation 
functions (S^^S^) (/i = x, z) and the quantities which can be obtained from {Sj J '(t)Si +n ) 
by partial or complete Fourier transformation with respect to space and time, namely the 
spatial Fourier transform, which is sometimes also called intermediate structure factor, 

P(q,t) =J2eM-iqn)(Sttt)S? +n ), (1.2) 

n 

the temporal Fourier transform 

/CO 
dteM^t)(S?(t)S? +n ), (1.3) 
-oo 

and finally the dynamic structure factor 

/oo 
dteM-i(Qn-cot))(St(t)St +n ), (1.4) 
-oo 

which is related to the inelastic neutron scattering cross sectionS. We have determined 



the above quantities by complete diagonalization of the Hamiltonian ( |1 . 1| ) for system sizes 
iV < 16 and anisotropies A = 0,cos0.37r, and 1, at various temperatures. Given the 
eigenvectors |z/) of Hxxz and the corresponding eigenvalues E u , the dynamic structure 
factor may be computed as 

S"(q,u) = ^J2 e M-PE„)5(u -E x + ^)|(z/|^|A)| 2 , (1.5) 

v, A 

where 

Z = 5>xp(-/^) (1.6) 



1 N 



SZ = ^=Y,exp(-iql)S?. (1.7) 



is the partition function, and 



N l=1 

The plan of the paper is as follows. In Sec. [Tl] we give a survey of previous (analytic) 
results, for the XX model ( [11 A| ) and for the XX Z model ( [11 B| ), we also recall the phe- 
nomeno logical pictures of ballistic and diffusive dynamics ( |11 C| ). Sec. [Ill A] addresses the 
question, which matrix elements give the largest contributions to the dynamic structure fac- 
tor as defined in Eq. ( |1.5| ). We have studied this question in some detail for the XXX- model 
(A = 1) and, for purposes of comparison, also for the frustrated Heisenberg model and the 



Haldane-Shastry model. In Sec. |III B| we present our results for the dynamic structure fac- 
tors. Sec. [Ill C| contains remarks on finite-size effects. Sec. |IV| is devoted to time-dependent 
correlation functions and Sec. [V| contains a summary. A detailed comparison of our numer- 
ical results (including also effects of a nonzero magnetic field) to recent experimental data 
from spin-chain compounds like KCUF3 or CuGeOs is planned for a later publication. 



II. PREVIOUS RESULTS 



A. The XX model 



For the XX model (A = 0), which can be mapped to a system of noninteracting lat- 
tice fermionsi'i by the Jordan- Wigner transformation, time-dependent correlation functions 
have been calculated^ at zero and nonzero temperatures. In that case the spin correla- 
tion function (Sf(t)Sj) is a simple fermion density correlation function, and the function 
(Sf(t)Sj) can be reduced to a determinant whose size increases linearly with i + j. 

Niemeijerd and Katsura et al.i derived explicit expressions for the correlation (Sf(t)Sf +n ) 



of the XX model and the associated quantities ( |1.2[ - |1.4|) which we will recall below. 



For arbitrary T and infinite N the structure factor is 



e(|4Jsin(f)|-M) 



2[(4J) 2 sin 2 (f) -u 2 ]? 
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exp(^; 



cosh 



& + / gjoos(|)[l 



(4J)2sin^(f) 



IT J 
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— [3 J cos( 



[1- 



I ^" ^"\2>^ (4/) 2 sin 2 (§) 



11/2 



(2.1) 



At the upper boundary of the frequency spectrum S z (q,u) displays an inverse square- 
root divergence for all temperatures. This divergence translates into an asymptotic 
t -1 / 2 exp(i4Jt sin(q , /2)) behavior for I z \q,t) and the existence of a high-frequency cutoff 
also leads to the long- time asymptotic behavior (S[(t)S z +n ) ~ t~ l . At T = 0, the dy- 
namic structure factor S z (q,u) as given in ( |2.1| ) has also a lower spectral boundary, where 
it displays a discontinuity. The longitudinal spin pair correlations at infinite T are 

1 



(2.2) 



an 



J 2 (g,t) = ij (4Jtsin|). 



(2.3) 



The transverse correlation functions (Sf(t)Sf +n ) of the XX chain are complicated many- 
particle correlations of the Jordan- Wigner fermions and fewer explicit results are available 
than for the longitudinal correlations. In fact, the only simple closed- form result is0~@ 

,2 



UJ 



for T = oo. Its et al.El determined (Sf(t)Sf +n ) for t — > oo, 
moderate magnetic field. 



n 



(2.4) 

oo and n/t = const in a 



B. The XXZ model 



At T = the long-distance and long-time properties of the XXZ chain may be de- 
scribed by a Fermi field theory0'0. This leads to power-law singularities with A-dependent 
exponents for several quantities. Specifically we have 
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W)sf + j~(-i)> 2 -^ 2 r^ (2.5) 



with 



*3b = *77 1 = l-^ (2-6) 

7T 



where 



A = cos7, (2.7) 



and 



ii J sin 7 . 
«= -. 2.8 

7 

This leadsS to corresponding low-frequency singularities in S^(q,u) and $q(cg>). The struc- 
ture of these singularities has been the subject of a recent studyEl. 

The continuum of (q, u) values for which S z (q, u, T = 0) is nonzero in the XX model is 
a special case of a more general continuum in the XX Z case. The spectral boundaries of 
this two-spinon continuum are given by 

7rJshi7 

Ei{q) = sing, (2.9a) 

7 

. . 2-7rJsin7 . q . ^, . 

e 2 (g) = -smj. (2.9b) 

7 2 

The nature of the two-spinon states defining the continuum will be discussed in Sec. [Ill A| . 
It was arguedS that at T = the dominant contributions to both S x (q,u) and S z (q, uS) 
come from this continuum and a related one with the same lower boundary and the upper 
boundary 

e 2 (q) = e 2 (jr-q). (2.10) 

Approximate analytic expressions for the T = dynamic structure factors were conjectured, 
which take into account known sum rules as well as exact results for the case A = and 



6 



asymptotic results Q2.5Q for small q and u. For explicit formulae and further discussion we 
refer the reader to Ref. |19[ Very recent lyEHO] the exact two-spinon contribution to S IJ '(q : oj) 
of the Heisenberg model has been be calculated analytically at T = and it was found that 
it accounts for more than 80 percent of the total intensity. 

Schulzi used field-theoretical techniques to study a spin-S" XX Z antiferromagnetic chain 
at finite (not too high) T. He derived approximate analytic expressions for S^^q, u) valid for 
small values of u and \q — 7r|, which showed good agreement to recent T-dependent inelastic 
neutron scattering datall from the 5 = 1/2 Heisenberg antiferromagnet KCUF3. For the 
case of interest, the results for 5^(5, o>) derived by Schulz may be written as 

S"(q,u) = ^( 7 )K + l)Im(sin(^)r 2 (l - \) 



2naT 



{ y 2 ' y 2 

p{x + ) P {x-) - ^ j . (2.11) 

Here 

= (exp(^) - I)" 1 , p(x) = T ^-™ ) ix) , x± = U± ll q T ~ n) (2.12) 



and v and rj^ are defined in fl2.8|) and (2.6). a is a short-distance cutoff (comparable to the 



lattice constant), and ^(7) is an overall factor. The <5 MZ term serves to cancel a spurious 
divergence which occurs for /i = z as 7 — > |. At T = 0, (|2.11 ) yields 



S»(q,u) ~6(cu- \vq\){u 2 - {vq) 2 ) 2 1 . (2.13) 

This expression correctly reproduces the T = singularitiesEl'Ellli at the lower continuum 
boundary £\{q) ( |2.9a| ). Note, however, that the experimentally confirmed!! upper continuum 
boundary s 2 (q) (|2.9b| ) does not show up in ( p.l3j ). Consequently the frequency integral over 



S IJ '(q,uj) diverges for arbitrary q and 7. 

Schulz' field theoretical approximation ( |2.11| ) reads in case of the XX-model and q = 71 



S z (q = vr, co, T) oc l — R (2.14) 

^ ' 1 (l + exp(^)) 2 1 } 



2 

which is the exact result Eq. (|2.1|) up to the factor 



e ( 4J -H), (2.15) 

[(4J) 2 -cu 2 ]5 

That is, Schulz' formula describes the T dependence of the exact XX result correctly, 
but fails to reproduce the cutoff and the square root singularity in u> at the cutoff. 

C. Phenomenological pictures: Ballistic and diffusive behavior 

Spin diffusion has been a popularii concept used in the discussion of experimental results 
for decades. According to the spin diffusion hypothesis one expects 

(St(t)SU - (S?S? +n ) ~ n exp(-— ) (2.16) 

for long times and long distances, where D is the diffusion constant. The intermediate 
scattering function correspondingly is 

P(q, t) - I"(q, 0) ~ exp(-DqH) (2.17) 

for small q and long times. 

A quantity sensitive to spin diffusion, the "spatial variance", was introduced by Bohm 
and Leschkell. It is defined as follows: 

= 4E(W*)^fiJ - WS&J)n 2 . (2.18) 

n 

Specifically, cr 2 (t) allows to distinguish between diffusive and ballistic behavior of the spin 
correlations. In the diffusive case, the intermediate structure factor ( |2.17|) yields 



- ^(n^)-n?,o)) 



= ioJ(t) ~ 2Dt (2.19) 

q=0 



On the other hand, ballistic propagation leads to 

o*(t) ~ t 2 (2.20) 

with the mean square velocity as a factor of proportionality. 
For the XX chain and T = oo one gets 
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al(t) = (2.21) 

and from (|2.3|) one obtains, by (|2.19|) 

a 2 (t) = 2JH 2 , (2.22) 

meaning that in that special case the conserved z spin component propagates ballistically, 
whereas the non-conserved x component decays locally. At finite T one obtains 

Ree^(t) oc t 2 , lma 2 (t) oc t (2.23) 

for A = and arbitrary t. It should be noted that for short times, ( |2.23|) is the leading-order 
behavior to be expected for both u 2 {t) and cr 2 (t) (and for zero as well as nonzero A) due to 
the general symmetry {Sj 1 {t)Si +n ) = (5f(— £)S^ n )*. From the asymptotic behavior found 
by Its et al.@ a 2 x {t) is expected to show exponentially damped oscillatory behavior for long 
times. 

Numerous theoretical studies have attacked the question whether spin diffusion describes 
the high-temperature dynamics of quantum spin chains; we will quote only a few selected 
references. In an important early numerical study, Sur and Lowe@ computed the XXX spin 
autocorrelation function at T = oo for N < 11. A comprehensive T = oo moment study0 
was extended a few years lateiS It was conjecturedil that spin diffusion is absent in the 
XXZ chain and numerical evidence to that effect was reported recently^. Other recent 
publications employ different approximate methods and arrive at different conclusions. In 
Ref. [3(] spin diffusion is not found, whereas in Refs. |3~T1 , |32"| spin diffusion is found. We have 



studied the longest chains to date, and for the time range covered in our study, the T = oo 
results presented in Sec. |V| do not show the signatures of spin diffusion. 



III. DYNAMIC STRUCTURE FACTOR 
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A. Excitation continua and classes of states 



Here we shall first remind the reader of the main conclusions at T = and then describe 
the corresponding results for T > 0. 

For the T = spin dynamics a very simple picture has emerged : only transitions which 
excite not more than two spinons are of importance. This statement will be made more pre- 
cise below. In the following we have to use the Bethe-Ansatzi! terminology (see appendix). 
Excitations with a small number m of flipped spins with respect to the ferromagnetic ground 
state have spin wave character and are called magnon excitations. They are relevant for the 
ferromagnetic model and for very strong magnetic fields in the antiferromagnetic model. For 
excitations of the antiferromagnetic ground state with m ~ N/2 and large N a description 
in terms of magnons is inappropriate. These states are fully characterized by the posi- 
tions of holes in their A - sequence. Excitations of this kind have first been studied by des 
Cloizeaux and PearsonB The first complete determination of the two parameter continuum 
corresponding to two independent holes in the A sequence of S = 1 states was presented 
by Yamada0 and Miiller et al.0. A systematic description was later given by Faddeev and 
Takhtajan0'il. Some important examples were worked out for the anisotropic model by 
Woynarovichil. These authors found that the elementary excitations of the antiferromag- 
netic Heisenberg model have spin 1/2 and that their number is always even. In recent years 
they were named spinonsS. The momentum of a spinon has the range < k < tt and its 
energy is 

e(Jfe) = j7rsin(A;)^M (3.1) 

7 

The momentum and energy of an eigenstate of Hxxz relative to the ground state are the 
sum of the momenta and energies of the individual spinons building up the state. Two- 
spinon eigenstates thus yield the continuum described by fl2.9a| , [2.9fj|) . Every eigenstate of 
Hxxz belonging to the set Maf (see appendix) has a fixed number of spinons. States 
with spin S have spinon number N sp = 2S + 2k, k = 0, 1, ... For finite S the class C (see 
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appendix ) of states turns out to be identical with the set of states having the smallest 
spinon number for fixed spin S and maximal S z = S. The ground state in the subspace with 
spin 5* is always a class C state with spinon number N sp = 2S. For T = and magnetic 
field h > one observes for matrix elements of local spin operators the selection rule, that 
the spinon number changes by and ±2 only. This rule is approximately valid for the 
Heisenberg model and exact for the Haldane-Shastry model. But by far not all transition 
matrix elements obeying this rule are large. If two states differ in spinon number by or 2, 
their energy difference may be much larger than €2{q) and in that case the matrix element 
is always small. 



This has been observed in Ref. [36] and confirmed by us for larger systems. The dominant 
ground-state matrix elements for magnetic field h > with spin S come from two parameter 
continua of states which are related to the ground state by the simplest hole excitations in A 
space (the classes SWC1 and SWC2 defined in Ref. ^ for longitudinal correlations), whereas 
the dimensions of the continua defined by AN sp = 0,2 are 2S, (2S + 2) respectively. That 
means that only a small subset of all transitions with AN sp = 0,2 contributes substantially. 

In the following we shall describe the nature of the excitations relevant for T > 0. We 
find that S z (q,u,T) is for all T almost completely confined to 

M < e 2 (<7), (3.2) 

which is the 2-spinon boundary ( |2.9b|) . The order of magnitude of suppression of the states 
lying outside of these limits can be read off from FigfJ for low T and from Fig.|| for high T. 
We reached this conclusion by 

1. a detailed examination of a set of selected states like those from the S = two-spinon 
continuum degenerate in energy with the S = 1 two-spinon continuum and also higher 
class-C states with S — 1, 2, 3 . In all these cases we found excitations with AN sp < 2 and 
u < e 2(<?) to be dominant. 

2. making a search in the set of all matrix elements (approximately 22.3 million ) to 
detect sizable contributions with exceptionally large u values. Not a single one was found. 

11 



The selection rules found for T = imply that only class-C states are important : 
as ground states for given external magnetic field and as excited states. Therefore it is 
interesting to assess the role of class-C states for T > too. 

We computed for the XXX-model for N = 12 and 16 the energies of all class-C states 
using the Bethe ansatz. To describe the results we introduce an additional label a in the 
symbol for a state \u,a) where a indicates the class to which the state belongs: a class-C 
state and all elements of its spin multiplet have a = 1, a bound state (see appendix) has 
a = 0. We define S aia2 (q, u>, T) by a formula equivalent to Eq. ( [L.5D with \u), |A) replaced 
by \u, ax), |A, CK2), and fixed ot\, 012. By integration over u we obtain the four quantities 
5*00 ; Sqi, Siq, S\\. The result is shown in Table | for T = 00. It is evident that for T 7^ and 
growing N bound states become increasingly important. 

In Table |J we show further characteristic quantities derived from the numerical results 
for S z {q,uS) at T = 00 for 7 = and 7 = 0.37T. For each q we list the upper continuum 
boundary £2(9) Q2.9b|) and the fraction 



XouM) = ' , Qz , v 3.3 

of the spectral weight situated outside the continuum. Even for small q, where £2(9) is small, 
x out (q) is not large, and it becomes rapidly smaller as q grows. As q — > ir the fraction x out (q) 
becomes negligible for all practical purposes. 

The excitation continua governing the longitudinal and transversal dynamical structure 
factors in the anisotropic Heisenberg model will be described in detail elsewhere. 

We conclude this section with results for non- Bet he- Ansatz models. The exactly solvable 
Haldane-ShastryliiS model given by the Hamiltonian 



H HS = 2J £ d(m - n)' 2 S m ■ S n (3.4) 

with 

N Ten 

din) = — sin—— (3.5) 
71 N v ' 
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describes an ideal spinon gas0. The dynamical spin correlation functions are known exactly 
for T = (see Refs. but there are still no complete results for T 7^ (see Ref. pE5| ). 



We determined for this model all transition matrix elements as given in Eq. ( |1.5| ) for 16 
spins. Only u which are multiples of n 2 J/N 2 occur. For small chain lengths S z (q, 00) is far 
from its continuum form because the energy eigenvalues are highly degenerate. We plot in 
Fig.||] the sum of all squared matrix elements for the discrete set of equally spaced u. For 
iV = 16 we find that the upper limit for uo given by Haldane and Zirnbauer@ for ground 
state excitations 

^max = ^"<?( 27r - q) (3-6) 

is strictly valid for all matrix elements contributing to S z (q,u,T). This completes observa- 
tions made by TalstraEl who studied a large number of examples. 

Finally we discuss results for another non-Bethe-Ansatz model. To explain the properties 
of quasi-one-dimensional compounds the description in terms of a spin model including inter- 
actions between both nearest and next-nearest-neighbors may be necessary. The frustrated 
model 

N 

Hfr W = 2J^ {Si ■ S i+ i + aSi ■ 6*4-2}, (3.7) 
i=i 

46| 



has been used in attempts to describe experimental data for CuGeOa with a = 0.24c9 and 
Sr 14 Cu240 4 i with a = 0.5 (see Ref. |47|). It is therefore interesting to find out whether the 
spin dynamics is significantly changed when additional interactions are introduced. In Fig.|] 
we compare S z (q = ir, u, T = 00) for two a values. We see that the spectral weight vanishes 
outside the two-spinon boundary of the isotropic Heisenberg model for both models. 

B. Results for S ll {q,u) 

In this section we describe the results for the dynamic structure factors S tl (q, uj) for Hxxz 
(|1.1[). For any finite N, S^(q,u) (as given in Eq.(|1.5|)) consists of a finite number of delta 
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peaks. Of course the number of peaks grows exponentially with N, but only a rather small 
fraction of them contributes significantly. The number of these relevant peaks varies strongly 
with N, as do their precise locations and strengths. (The regular structure visible in Fig.^Jfor 
the Haldane-Shastry model is a striking exception to this rule.) In order to avoid confusion 
by these finite-size effects we present our results for S^lq, u) as histograms, displaying for 
each frequency bin of constant size Au; the total spectral weight of the S^{q,uj) peaks in the 
bin, divided by Au. The size of the frequency bins chosen represents a compromise: too 
narrow bins overemphasize finite-size effects, too wide bins erase all structure. 



1. The XX model 

For the XX model, S z (q, uj) is known in closed form (see |2.1| ), but S x (q, uj) is only known 
numerically, except at infinite T (see [2.4] ). The shape of S z (q,uj) (see Fig.|5|) is dominated 
for low T by the spectral boundaries at e\ and e<i and in general by the inverse-square-root 
singularities uj = ±£2(9) ( |2.9b| ) . For finite T the cu-dependence of S z (q, to) thus shows a 
typical "trough" shape (see Fig.|). This trough shape is a feature which is present also in 
the XX Z model results to be discussed below. 

The function S x (q,uj) as shown in Fig.^ involves in the Jordan- Wigner fermion picture 
the excitation of arbitrarily many particle-hole pairs and thus its frequency range is not 
restricted. This is reflected in the g-independent Gaussian shape of S x (q,uj) at T = 00: 
spectral weight is present at all frequencies, but it is very small for high frequencies. The 
"hill" shape exemplified by the Gaussian again is a generic feature which persists also for 

7 7^ 0.57T. 

At 7 = 0.57T we find the typical hill structure in S x (q,u) for all q and for all T > 2 J. At 
low T the bulk of the spectral weight is contained in the region confined by the continuum 
boundaries £1(5), £2(5), and £2(9) conjectured^ for T = 0. 
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2. The XXZ model, low T 



Here we discuss our low-temperature results, starting with the longitudinal structure 
factor. As discussed in Sec. [Ill A , the low-temperature N = 16 data for 7 = show that the 
spectral weight in S z (q,u) is almost completely confined to the continuum between £i{q) 
and £2(9)- 

Recently it was showniS that the two-spinon contribution to S x,z (q,u) vanishes in a 
square-root cusp at the upper boundary. 

The contributions to u) from outside the two- spinon continuum become visible on 
a logarithmic scale. In Fig.§ we have plotted S z (n/2, u) for A = cos(0.37r) and iV = 16. The 
two frequency bins containing almost all of the spectral weight lie within the boundaries of 
the two-spinon continuum; beyond £2(9)) S z (n/2,u) decays rapidly. 

At low T the continuum boundaries are generally quite well visible in the numerical 
results for both 7 = and 0.37T. (Compare Fig.^ for 7 = 0.37T and Fig.| for 7 = 0; q = 37r/4 
in both cases.) At 7 = and small q, there is some spectral weight rather far above the 
upper continuum boundary £2(9)- The maximum of S z (q, u) is close to the lower continuum 
boundary for those q values where the continuum boundaries are well separated on the scale 
given by the frequency bin width. For 7 = that is consistent with the conjecture of Ref. 



The low-temperature results for the transverse structure factor S x (q,u) at 7 = 0.37T 
(compare Fig.|] for q = 37r/4) are similar to those at 7 = 0.57T. The continuum boundaries 
are again quite well recognizable, with only very little spectral weight above the upper 
continuum boundary. The maximum of S x (q, uj) is always situated close to the lower spectral 
boundary and its height grows with q. 



We have carried out a detailed comparison between Schulz' field-theoretical result ( 2.11|) 
and our numerical N = 16 results for S^{q, uS) at q = 7n/8 and 7r. As expected, Q2.ll ) works 
best in the low-energy region at low temperatures. For T — > the singularities show up at 
the lower boundary £i(q) of the continuum. There are no sharp high-frequency cutoffs nor 
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singularities, but apart from that, the shapes of the two curves are similar. As an example 
we show in Fig.0 S z (7n/8, u) for the XXX model at T = 1.2J (see Ref. (4|). There the 
low-frequency data are described well by ( [2. lip , and the high-frequency cutoff is replaced 
by a gradual but quite rapid decay 



As the temperature is raised to T = 2.7 J (as seen in Fig.|TT|) Eq. ( [2. lip extends further 



beyond the range of the numerical data for both positive and negative frequencies. The 
maximum of (|2.11 ) has moved towards higher frequencies as compared to Fig.|lO| and as T 



is raised further the maximum moves to frequencies where the numerical data essentially 
vanish. 

The result (|2.11|) of Schulz was successfully^! used to fit the T-dependence of inelastic 



neutron scattering results on the Heisenberg antiferromagnetic chain compound KCuF 3 in 
the vicinity of q = tt. We wish to point out that the agreement with the experimental data 
found in that study does not contradict the differences to our numerical results described 
above. Firstly, the experiments were all carried out in the range T < J where ( [2. lip is 



very good at low frequencies and not unreasonable at high frequencies. Secondly, ( |2.11p was 
only used to describe the experimental data in the region q ph tt at frequencies far below 
the upper continuum boundary e 2 (q). High-frequency data (see, for example, Fig. 17 of Ref. 
|23"1) were interpreted in terms of the T = conjecture of Ref. [19|, and they did show a clear 
spectral cutoff consistent with that conjecture (and thus also with our numerical data). 

3. The XX Z model, high T 

The temperature dependence of the dynamic structure factor displays clear general 
trends: For small q, 5* M (g, to) is largest at high T, whereas at larger q, S^g, u;) decreases 
with T. For A ^ 0, the high-temperature dynamic structure factor S^(q, u) always has its 
maximum at small q. Given that S z (q,u,T) is almost completely confined to \u\ < £2(9) 
this is a consequence of the sum rule 

J dcuS z (q, w, T = 00) = tt/2. (3.8) 
16 



The general shape of the high-temperature structure factor S IJ '(q, uj) (for iV = 16) changes 
smoothly as one goes from 7 = 0.57T to 7 = 0: at 7 = 0.37T S x (q, uj) shows the "hill" shape 
familiar from 7 = 0.57T, and S z (q,uj) shows the "trough" shape. At 7 = 0, S z (q, uj) is of 
"hill" type for q < n/2 and of "trough" type for q > 71/ 2. 

Clearly, the number of terms contributing to S^^q, u) ( |1.5| ) at T 7^ is much larger 
than at T = 0, and as T — > 00, the contributions of different states \u) in (|1.5| ) are not 
discriminated by the thermal weight factor exp(— (3E V ) any more. Nevertheless Figs.[|-f| show 
that S^(q, uj, T) is confined to a narrow interval in uj, as discussed in Sec. [Ill A[ In Fig.f?] we 
show S z (^-,uj) at 7 = 0.3vr for a range of temperatures. The continuum boundaries (|2.9a|) 
and ( fT9Ep are e x = 1.906 J and e 2 = 4.982 J. For low T, S z (^-,uj) is basically confined to the 
interval between E\ and e 2 and shows sharp structures (probably due to finite-size effects) 
and a maximum close to E\. For higher T these structures are washed out and the probability 
of negative-frequency contributions to S z {^f,uj) grows in accordance with detailed balance 
S^(q, —uj, T) = S^(q,uj,T) exp(—/3uj) . No significant spectral weight appears for \u\ > e 2 
as T grows. 

Fig.0 shows S z (^,uj) at T = 00, for 7 = 0,0.3vr (for N = 16) and 7 = 0.5vr for 
iV = 00. The maxima of S z {^f,oj) for 7 = and 0.37T coincide with the singularities for 
the case, 7 = 0.57T. For frequencies beyond the maxima (where S z {^f,uj) vanishes strictly 
for 7 = 0.57r) we see a sharp decline for 7 = and 0.37T with small tails extending outward, 
roughly to £2(^1 7) ; where S z {^f, uj) becomes negligible on a linear scale. (On a logarithmic 
scale, rapidly decaying contributions for large |a>| are visible, similar to those in Fig.|]; see 
also Fig.|15|.) The observations from Fig.|12| illustrate the data presented in Table |I[ 

The transversal dynamical structure factor S x (q,uj,T = 00) evolves for varying 
anisotropy (from 7 = to vr/2) from a strongly g-dependent function with spectral bound- 
aries ±62(9) (see ( [2.9 b| )) to the g-independent Gaussian Eq. ( |2.4| ) (see Figs. |I],|T3"| for details). 
There is no physical explanation for this change of shape yet. 
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C. Some remarks on finite-size effects 



In contrast to the case of static structure factors, where systematic finite-size analyses 
can be performed at both zero and finite T, the presence of an additional variable makes the 
situation more complicated for dynamic structure factors. As stated at the beginning of this 
section, the number, the weights and the positions of the 5 functions contributing to S IJ '(q, to) 
( |1.5| ) vary with the system size N in a complicated manner. In the histogram representation 
which we use this means that as N varies, spectral weight is exchanged between neighboring 
bins in a seemingly random manner. 

For the rather wide bins used here and for high T, where many states \v) in ( |1.5|) 
contribute the weight exchanged between two bins is very small compared to the total 
weight in each bin and thus finite-size effects are small. At low T, where a small number 
of 5 peaks dominates S^lq, to) we do indeed observe larger finite-size effects. This is seen in 
Fig.|14| (main plot), where we show S z (tt,uj) at T = 0.2J, 7 = 0.37T for N = 12, 14, and 16. 



The same quantity, at T = 2.2J is shown in the inset of Fig.|l4] demonstrating that with 
the general smoothing effect of nonzero T, the finite-size effects in S^(q, u) also become 
smaller rapidly. The rapidly decaying spectral tails in S^(q,u!) which extend well beyond 
the two-spinon continua show only very little size-dependence at T = 00, as shown in Fig.|I~5| 
for q = 7r/2 at 7 = and 0.37T. 

IV. TIME-DEPENDENT CORRELATION FUNCTIONS 



In Fig.[16] we show the real part of (Sf(t)Sf +n ) for n = 0,1,..., 8 at T/J = 20, for 
N = 12, 14, and 16, in the XX Z chain with A = cos(0.37r). As in all plots of this type the 
results for different N are identical at short times and differ from each other at longer times 
in a systematic way, with the data for the smallest N deviating first. (In the XX case this 
behavior can be observed over a wide range of system sizes.) Thus the results are certainly 
valid for N = 00 over the time range where N = 12, 14, and 16 yield identical correlations. 
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The onset of finite-size effects occurs earlier for larger n. This is easily interpreted by 
observing that the finite-size effects are due to correlations connecting the sites I and I + n 
"the wrong way round", i.e. over a distance N — n instead of n. The (roughly) linear growth 
with n of the time range affected by finite-size effects is consistent with ballistic propagation 
of the correlations. In an analogous way one may understand the linear growth with n of 
the short-time region of practically vanishing correlations: the information from site / has 
not yet arrived at site / + n. 

At T = oo, we compared the finite- N data to exact infinite- N resultsi! obtained by 
moment-expansion methods. These methods employ the power-series expansion coefficients 
of (Sj l +n (t)Si') (which is an even real function of t at T = oo) up to a certain order to obtain 
exact upper and lower bounds to (S'/l | _ Tl (t)S'f ) (or to I^(q,t)) for the infinite system. Over a 
certain time range the upper and lower bounds coincide (for all practical purposes) and then 
they start to deviate from each other rapidly. Bohm and Leschke§! obtained the expansion 
coefficients to order t 30 for the XXX-chain and the resulting bounds coincide with each 
other (and thus with the exact N = oo result) for Jt < 2.8. In Fig.|H] we show these bounds 
together with our results for N = 14 and N = 16. 

In Fig.[T3| we display the z spin autocorrelation function (Sf(t)Si) for A = cos(0.37r) 
and four representative T values. Two general tendencies can be observed in these data as 
T is lowered: firstly, the oscillations in the correlation function become more pronounced, 
and secondly, finite-size effects start to show up earlier at lower T, due to the small (size- 
dependent) density of states at low frequencies. The dominant time scale of the oscillations, 
however, does not change with T. 

In Fig.|19] (main plot) we show the real part of the transverse autocorrelation function 
(Si(t)Si) of the XX chain for three (high) temperature values, in comparison to the exact 
Gaussian infinite-T (and infinite- N) result. At short times, all results coincide. The T/J = 2 
results branch off first, the T/J = 20 results follow later. The T = oo results for iV = 
16 coincide^ with the Gaussian over the time range shown here. In the non-Gaussian 
regime (for finite T), the data for N=16 and 14 stay together for a while, before finite-size 
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effects become important. The short time interval between the Gaussian and finite-size 
dominated regimes displays the exponential decay predictedlll and verifiedEl recently, with 
the expected value of the decay rate. The inset of Fig.[19] demonstrates the changes that occur 
in the infinite-T x autocorrelation when the interaction A is switched on. A finite A both 
accelerates the initial decay of (Sf(t)Sf) (an effect well known from short-time expansion 
studies at T = oo0'i!) and decelerates the long-time decay, which becomes much slower than 
the exponential decay present at A = 0. This observation is evidence for the exceptional 
character of the high-temperature x autocorrelation function at A = 0. The limited time 
range covered by our data prevents a quantitative characterization of the changes induced 
in the long-time behavior by a nonzero A§. 

In contrast to the rather pronounced slowing-down of the long-time decay of (Sf(t)Sf) 
with growing A at high T (see Fig.|l^ (inset)), the long-time decay of (S[(t)S[) varies only 
weakly with A. This is demonstrated in Fig.^0] (main plot) for T = oo: the zeros of 
(Si(t)Sf) at A = (compare ( |2.2j )) gradually "fill up" as A grows, but the overall decay 
stays roughly the same. The curves for Ke(Sf(t)Sf) at T/J = 20 are very similar. At 
T/J = 0.2 all autocorrelations show much more structure than at high T, but neither the x 
autocorrelation (not shown) nor the z autocorrelation (Fig.^, inset) show strong variations 
of the decay as A changes. 

We now discuss the "spatial variance" cr 2 {t) as defined in ( [2.18|) . For A = 0, cr 2 (t) is 



given by ( f2.23|) , and that is precisely what the numerical results show, apart from finite-size 



effects. For given N and T, Rea 2 (t)/t 2 and lma 2 (t)/t are constant up to some characteristic 
time t c which is approximately proportional to N and grows slightly with T; at high T we 
have t c m N/4 J. For t > t c finite-size effects start to build up. Assuming that t c is the time 
it takes the fastest excitations to travel the distance N/2, one obtains a velocity v — 2J 
for these excitations. This is the maximum velocity of the Jordan- Wigner fermions and 
also the spinon velocity ( |2 . 8| ) for 7 = ir/2. The root-mean-square velocity [Kea 2 (t) /t 2 ] 1 ^ 2 
(compare ( |2.20[ )) is smaller (by a factor a/2, see ( 2.22j) ) than the maximum velocity at high 



temperatures. With decreasing T, Kea 2 (t)/t becomes smaller, whereas lma 2 (t)/t grows, as 
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may be seen by combining (|2.19 ), and a high-temperature expansion to obtain 



2 t 2 
and 

IwJ(t) 



' 2 ^Kt~\ (4.1a) 



ocT- 1 . (4.1b) 



t 

At A = cos(0.37r), Recruit) /t 2 and lma 2 (t)/t are not constant, but decrease slightly from 
their t = values, until finite-size effects set in. This decrease becomes stronger at A = 1. 
The dependence of the t = values of Kea 2 (t)/t 2 and \ma 2 {t)/t on T is roughly similar 
for all three A values considered here (see ( f4. Ia| , f4. ) However, \im t ^ Real{t)/t 2 is a non- 
monotone function of T at A ^ 0, with a maximum at intermediate T; at high T it becomes 
roughly T- independent. 

The spatial variance cr 2 {t) of the x correlation shows qualitatively the same behavior for 
all T and A values studied here. 'Rea 2 {t)/t 2 and lvaa 2 {t)/t both decrease from their t = 
values, with the decay becoming more rapid as one goes from 7 = 0to7 = 7r/2. The 
t = values of these quantities themselves drop rather rapidly as 7 grows. For 7 = 0.37T the 
T-dependence of lim^o R- ecr xWA 2 i s similar to the 7 = case, starting with a small value 
at low T and saturating at a higher value for high T, whereas at 7 = 7r/2 lim^o R- ecr x(^)/^ 2 
starts out with a small value at low T and then rapidly drops to 0. To summarize, the data 
for <J 2 (t) do not provide evidence for diffusive behavior within the time range accessible to 
our study. The behavior of <J 2 (t) does not deviate dramatically from the form ( |2.23| ) valid 
for A = at arbitrary times and for arbitrary A at short times. 

We finally discuss the intermediate structure factors I x (q,t) and I z (q,t) as defined in 
( |1 . 2| ) . In the XX model we find at T / J = 20 a g-independent decay of I x (q, t), very similar 
to the T = 00 Gaussian. The deviations from the Gaussian decay which become visible at 
lower T were already discussed above in terms of the spin pair correlations. Similarly, the 
evaluation of I z (q, t) for the XX model yields no new insights as compared to the spin pair 
correlation or the dynamic structure factor. 
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For the general case of the XX Z model, however, both I x (ir,t) and I z (ir,t) show an 
exponential t dependence for sufficiently large t and T values!!. This is shown in Fig.|2T|, 
where we have plotted the absolute values of I x (7T,t) and I z (ir,t) for T = oo . 

In the time range Jt < 5 the T/J = 20 curves coincide with the T = oo curves. At 
T/J = 2, the exponential decay is much less marked and for even smaller T it is no longer 
perceptible. 

V. SUMMARY 

We have performed the first calculation of dynamic correlation functions for the antiferro- 
magnetic spin- 1/2 XX Z chain (|1.1|) at arbitrary temperature by complete diagonalization 
of systems with size N < 16 and anisotropics A = 0, cos(0.37r), and 1. We have calcu- 
lated the dynamic structure factors S IJ '(q, uj) (// = x,z) (|1.4j) and their Fourier transforms 
(S?(t)S? +n ) and I»(q,t) flL 



Our study extends the (T = 0, N < 10) work of Refs. ^§,|T^ to larger systems and to T > 
0. We find that even at nonzero T only a very small fraction of the many possible transitions 
yield appreciable contributions to S fl (q, uj) (see Sec. [Ill A|) . At all temperatures up to T = oo 



the dynamic structure factor S z (q,u) turns out to be of negligible size for greater than 
the upper boundary of the knowiffl excitation continua. This is demonstrated in Figsffl-M 
Table H shows the the very small fraction of spectral weight of S z (q, uj) outside the two- 
spinon continuum at T = oo. There are analogous results for the transversal structure factor 
S x (q, u) but this function becomes g-independent for A = 0. Similar conclusions about the 
spectral range of S^(q,uj) hold also for the frustrated Heisenberg model (|3.7|) (see Fig.|]), 
and for the Haldane-Shastry model (|3.4| ) (see FigfJ), where the T = oo dynamic structure 
factor for N = 16 vanishes strictly outside the two-spinon continuum. 

The general shape of S IJ '(q, uj) at high temperatures resembles that of the (analytic) 
T = oo results for 7 = 0.57T, where S x (q, uj) is a (g-independent) Gaussian "hill" and 
S z (q, uj) is defined on a restricted frequency range (which grows with q) with inverse square- 
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root divergences at the boundaries which lead to a distinct "trough" shape. These basic 
patterns are encountered again at 7 = 0.37T. At 7 = 0, S^(q, to) is a "hill" for q < n/2 and 
a "trough" for q > 7r/2. 

The low-temperature dynamic structure factors for small u and q ~ rr are described 
very well by the field-theory result ( |2.11D of SchulJl (see Fig.[Tt]) which contains no upper 



frequency cutoffs and can not be used at high temperatures (see Fig. [TTl) . 

Some results for space- and time-dependent correlations (Si(t)Si +n ) are discussed in 



Sec. [TV]. These correlations turn out to be iV-independent over a time interval growing 
with N. For 7 = 0.57T the long-time asymptotic behavior of (Sj J '(t)Si +n ) is known (see Sec. 
I|D, and under suitable circumstances that known long-time behavior may be detected in 
an N = 16 system (Fig.p!9|). However, reliable conclusions for the long-time asymptotic 
behavior (including the absence or presence of spin diffusion) at 7 < 0.57T cannot be drawn 
from iV < 16 data. For sufficiently high T and large q, the intermediate structure factor 
^(q, t) for 7 = and 0.37T displays an exponential time dependence (Fig.pTp. 
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APPENDIX A: BETHE ANSATZ NOTATION 

The eigenstates of the Hamiltonian ( |Lip expanded in the basis of states with fixed 
number r of down spins at positions x\, ...,x r have the coefficients 

g(x h ...,x r ) = 
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J^expi ( k ix x x + ... + k ir x r + ^Xj^wj • ( A1 ) 

The sum extends over all permutations i r of the integers 1, r. The pseudomomenta 
ki and phases ipji are solutions of the equations 

<PH A (cot y ' - cot |) 

cot ^ = i ? *2 A2 

2 l + A-(l-A)cot T cot| 

N/c m = 27rA m + <Pmi- (A3) 

A state is determined by the set of integers Ai, A r with < Aj < N. The corresponding 
energy eigenvalue is 

1 r 

E = -A(N -4r) + 2Y / cosk l (A4) 

2 i=i 

and the momentum is 

fc = |t^ (A5) 

JV m=l 

The ground state in the subspace with spin S (for A = 1) has 

Ai, X r = S + 1, S + 3, • • • , N - S - 1. (A6) 

The pseudomomenta fc, are real if Aj — Aj > 2. If in addition all Aj > the corresponding 
states are called class-C states!!. The remaining solutions are called bound states and have 
(with few exceptions) complex pseudomomenta forming (for large N) strings in the complex 
plane. The subspace Maf introduced i to describe spinon excitations contains all states 
for which y — V\ = finite for N — ► oo. z^i is the number of strings of length one. We 
have used Eqs.(|A2])-( |A5| ) to determine the complete set of class C solutions for A = 1 and 
iV < 16 for the analysis described in Sec. [Ill A. 
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TABLES 



N,g 


5ll 


Sio 


Sbi 


Soo 


16, 7T 


1.8990 TO" 1 


5.7144 TO" 2 


5.7144 T0~ 2 


6.9581 TO" 1 


12, TT 


3.5785 TO" 1 


8.8507 -lO" 2 


8.8507 T0~ 2 


4.6514 TO" 1 


16 - 


1.2673 TO" 1 


6.2628 TO" 2 


6.2628 TO" 2 


7.4802 TO" 1 


12 - 


2.4249 -lO" 1 


1.0750 TO" 1 


1.0750 TO" 1 


5.4251 TO" 1 



TABLE I. Spectral characteristics of the structure factor of the isotropic Heisenberg chain at 
T = oo. S a p is the fraction of S z (q,uj) corresponding to a transition from a class a state to a class 
(3 state, a = 1 denotes class C, a = denotes non-class C. 
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Q-- 


°2 VI) / ° 


XXX I \ 
^out ViJ 


°2 Vi)l " 


XX Z/ \ 

^out vy; 


1 


1.2258 


4 1539 -10~ 2 


1.0522 


4 5615 -10~ 2 


2 


2.4045 


2 4019 -10~ 2 


2 0639 


1 6524 -10~ 2 


3 


3.4908 


Q 5773 .10~ 3 


2.9964 


8 0768 -10 -3 


4 


4.4429 


4 1566 -10~ 3 


3.8137 


4 2199 -10 -3 


5 


5.2243 


1.7826 -lO" 3 


4.4844 


2.0314 -lO -3 


6 


5.8049 


7.5846 -10" 4 


4.9828 


8.9729 -10~ 4 


7 


6.1625 


3.4219 -10" 4 


5.2898 


4.1120 -10~ 4 


8 


6.2832 


2.4222 -10~ 4 


5.3934 


2.7909 -10~ 4 



TABLE II. Spectral characteristics of the structure factor of the XXX- and XXZ-Heisenberg 
chain (XXZ for 7 = 0.3) at T = 00: upper boundary of the two-spinon continuum and fraction of 
spectral weight outside the continuum. 
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FIGURES 
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FIG. 1. Dynamic structure factor of the N = 16 isotropic Heisenberg chain: logarithm of 
S z (q,uj) for T = oo and q = 7r/8,7r/2 and tt. 
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FIG. 2. The dynamic structure factor S z (q = 2 f E ,w) for the 16-site Haldane-Shastry model, at 
T = do. For each of the equidistant excitation frequencies we plot the sum of all squared matrix 
elements from (|1.5| ). Frequencies are multiples of 5 = 7r 2 /128 . The parabola-shaped spectral 
boundary as given by Haldane and Zirnbauer0 for T = (Eq. ^.6|) is strictly valid for all T. 
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OJ/J 

FIG. 3. Longitudinal dynamic structure factor S z (tt/2,uj) at T = of the N = 16 XX Z chain 
with A = cos(0.37r). The two-spinon continuum is marked by a horizontal bar at the 10 -3 level. 
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FIG. 4. The dynamic structure factor S z (tt,lo) for the 16-site frustrated Heisenberg model with 
a = 0.24 (solid line) and a = 0.5 (dashed line), at T = oo. 
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FIG. 5. Longitudinal dynamic structure factor S z (ir/2,oj,T) of the XX chain. (See Eq. 
Note the sharp cutoff at the lower spectral boundary near T = 0. 
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FIG. 6. Transversal dynamic structure factor S x (q = ir/2,u,T). The T = const cross section 
becomes an exact Gaussian for T = oo. 
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FIG. 7. Temperature dependence of the longitudinal dynamic structure factor S z (3tt/4, lo) of 
the N = 16 XXZ chain with A = cos(0.3vr). T varies from 0.2J (front) to 19.7J (back). The 
two-spinon continuum ranges from u = 1.907J to lo = 4.983J. Note that the x axis represents 
-co/ J for reasons of visibility. 
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FIG. 8. Same as Fig.|7|, for the transverse dynamic structure factor S x (3n/4, u). The bound- 
aries of the excitation continua are s\ = 1.907 J, £2 = 4. 983 J, and £2 = 2. 064 J. The x axis 
represents —to/ J for reasons of visibility. 
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FIG. 9. Temperature dependence of the dynamic structure factor S z (Sir / A, uj) of the N = 16 
isotropic Heisenberg chain. T varies from 0.2 J (front) to 19. 7 J (back). The two-spinon continuum 
ranges from 2. 22 J to 5. 803 J. The x axis represents —to/ J for reasons of visibility. 
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FIG. 10. Dynamic structure factor S z (7tt/8,uj) of the XXX-chain at T = 1.2 J. Shown are 
the numerical result from a N = 16 chain (solid line) and the field-theoretical result of Schulz, Eq. 



(2.11 ) (dashed line). 
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FIG. 11. Same as FigM for T = 2.7 J 
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FIG. 12. Longitudinal dynamic structure factor S z (3tt/4:, uj) of the N = 16 chain at 

T = oo, for A = 1 (dashed line) and A = cos(0.37r) (solid line). For comparison we also show 
the exact N = oo result for the XX chain, A = (thin line). The three pairs of triangles on 
the frequency axis mark the upper boundary ±e 2 (x) (H§) of the two 

-spinon continuum for 

A = 0, cos(0.37r), and 1 (in ascending order of \u\). 
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FIG. 13. Transversal dynamic structure factor: logarithm of S x (q, to) of the N = 16 XXZ-ch&va. 
with 7 = 0.37T and 7 = 0.57T at T = 00. a: g = | for 7 = 0.37T and the q-independent Gaussian for 
7 = 0.5-7T. b: 7 = 0.37T and q = 0,7r/2 and it. 
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FIG. 14. Finite-size effects in S z (ir,w) for 7 = 0.3vr at T = 0.2 J (main plot) and T = 2. 2 J 
(inset). Data for N = 12, 14, and 16 are shown as dot-dashed, dashed, and solid lines, respectively. 
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FIG. 15. Absence of finite-size effects at T = oo in the low-weight spectral tails. Main plot: 
S z (ir/2,uj) at 7 = 0.3-7T (solid line) and 7 = (dashed line). Inset: S x (tt/2,uj) at 7 = 0.3-7r(same 
range of variables as main plot). Data for N = 12 are shown as thin lines, those for N = 16 as 
heavy lines. 
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FIG. 16. Real part of the z spin pair correlation function (Sf(t)Sf +n ) for n = (top) through 
8 (bottom), at T/J = 20, for anisotropy A = cos(0.37r) and different system sizes, N > 2n, namely 
N = 16 (solid lines), N = 14 (dashed), and N = 12 (dot-dashed). (Curves are displaced in the 
vertical direction for graphical purposes.) 
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FIG. 17. The intermediate structure factor I z (q = tt, t, T = oo) for the XXX-chains of 14 and 



16 spins and the upper and lower boundaries (dashed lines) from Ref. 25. The curves for N = 14 
and N = 16 can not be distinguished on this scale. 
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FIG. 18. Real part (upper set of curves) and imaginary part (lower set of curves) of the z spin 
autocorrelation function (Sf(t)Sf) for temperatures T/J = 20, 2, 0.7, and 0.2 (top to bottom in each 
set of curves), for anisotropy A = cos(0.37r) and different system sizes, N = 16 (solid lines), N = 14 
(dashed), and TV = 12 (dot-dashed). (Curves are displaced in the vertical direction for graphical 
purposes; the y-axis scale applies to the uppermost curve; Re(Sf (0)Sf ) = 1, lm(Sf(0)Sf) = for 
arbitrary T.) 
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FIG. 19. Real part of the x spin autocorrelation function (Sf(t)Sf) of the XX chain (A = 0) 
for temperatures T/J = 2, 20, and oo, and system sizes N = 16 (solid lines) and N = 14 (dashed 
lines; no N = 14 data shown for T = oo). The N = 16 result for T = oo coincides with the exact 
result for N = oo following by Fourier transformation from (|2.4j ). The inset shows the N = 16 
results for T = oo, A = (dot-dashed), A = cos(0.37r) (solid line), and A = 1 (dashed line). 
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FIG. 20. Real part of the z spin autocorrelation function (Sf(t)Sf) of the XX Z chain for 
N = 16 and temperatures T/J = oo (main plot) and 0.2 (inset). 
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FIG. 21. The intermediate structure factor at T = oo. The solid curve repre- 
sents \I z (q = n,t, A = cos(0.37r))|, the dotted \I(q = n,t,A = 1)| and the dashed 
\I x (q = ir,t, A = cos(0.37r)|. For reasons of visibility the upper (lower) curve is multiplied by 
a factor of 10 (0.1). The straight lines are drawn through the fourth and fifth maxima. 
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